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I. INTRODUCTION 


According to the Schrodinger formulation of quantum mechanics, a total eigenfunction 
provides implicitly all relevant information about the behaviour of a physical system. Thus, 
if it is exactly solvable for a given potential, an obtained eigenfunction can be used to 
describe such a system completely. This has made the exact solutions of quantum systems 
be an important subject and also attract much attention in the development of quantum 
mechanics [1-25]. 

Up to now, there have been several efficient methodology developed to find the exact 
solutions of quantum systems within the framework of the non-relativistic and relativistic 



Very recently, Ma and Xu have proposed an exact (improved) quantization rule and 
applied it to calculate the energy levels of some exactly solvable quantum mechanical prob¬ 
lems. These include the finite square well, the Morse, the symmetric and the symmetric 
and asymmetric Rosen-Morse, the harmonic oscillator, the first and second Poschl-Teller 


potentials 


25| . By using this same rule, Dong and co-researchers 45j have obtained the 


analytical solutions of the Schrodinger equation for the deformed harmonic oscillator in one 
dimension, the Kratzer potential and pseudoharmonic oscillator in three dimensions. The 
energy levels of all the bound states are easily calculated from this quantization rule. 

The priority purpose of the present work is to obtain the obtain the bound state solution 
of Deng-Fan diatomic potential model via this quantization rule. The Deng-Fan diatomic 
potential model also known as the generalized Morse potential which was proposed some 


decades ago by Deng and Fan 


46[ |. This is done in an attempt to finding a more suitable 


diatomic potential to describe the vibrational spectrum, qualitatively similar to the Morse 
potential with the correct asymptotic behaviour as the inter-nuclear distance approaches 
zero. This potential model can be used to describe the motion of the nucleons in the 
mean field produced by the interactions between nuclei T7|. This potential has been used 
to describe diatomic molecular energy spectra and electromagnetic transitions and it is 
an ideal inter-nuclear potential in diatomic molecules with the same behavior for r ^ 0. 
Because of its importance in chemical physics, molecular spectroscopy, molecular physics 
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and related areas, the bound state solutions of the relativistic and non-relativistic wave 


equations have been studied by several authors 


0 


29 


46l. |47| |. The shape of this potential 


with respect to some diatomic molecules is shown in hgur^H 




FIG. 1: (color : online) Shape of Deng-Fan diatomic molecular potential for different diatomic molecules. 

The organization of this work is as follows. In the next section, we briefly introduced 
the exact quantization rule. In Section 3, we apply this quantization rule to obtain the 
bound state solutions of the Deng-Fan molecular potential. Section 4 present the numerical 
results and hnally and a brief conclusion. 


II. EXACT QUANTIZATION RULE 


In this section, the brief review of this quantization rule is given. The details can be 


found in Refs 


25 


26| . It is well known that in one dimension, Schrodinger equation 


dx^ 


i!{x) + 


2/i 

F 


[E - V{x)] ^p{x) = 0, 


( 1 ) 
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can be written in the following form 


+ (p{xY + k{xY = 0, with k{x) = — V{x)], 


( 2 ) 


where (j){x) = ijj\x)/'ijj{x) is the logarithmic derivative of the wave function 'ijj{x), the prime 
denotes the derivative with respect to the variable x, /i denotes the reduced mass of the 
two interacting particles, k{x) is the momentum and V{x) is a piecewise continuous real 
potential function of x. For the Schrodinger equation, the phase angle is the logarithmic 
derivative 0(a:). From equation ([2]), as x increases across a node of wave function 
(j){x) decreases to —oo, jumps to +cxo and then decreases again. 


In the recent years, Ma and Xu 25|, [2^ generalized this exact quantization rule to the 


3D radial Schrodinger equation with spherically symmetric potential by simply making the 
replacements x ^ r and V{x) —)■ Veff{r): 


I'rb rrh 

k{r)dr = Nir + / 0(r) 


dk{r) 

dr 


d(f){r) 

dr 


1 -1 


'=(>■) = 


Veff(r)]. (3) 


where and are two turning points determined by E = Veff{r). The N = n + 1 is 
the number of the nodes of (j){r) in the region Ene = 14 //(r) and is larger by one than 
the number n of the nodes of wave function -^(r). The hrst term Ntt is the contribution 
from the nodes of the logarithmic derivative of wave function, and the second is called the 
quantum correction. It is found that for all well-known exactly solvable quantum systems, 
this quantum correction is independent of the number of nodes of wave function. This 
means that it is enough to consider the ground state in calculating the quantum correction, 
i.e. 

Qc= [ Kir)^dr (4) 




In the recent years, this quantization rule has been used in man y physic al systems to obtain 


the exact solutions of many exactly solvable quantum systems 
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III. APPLICATION TO DENG-FAN MOLECULAR POTENTIAL 

In this section, we apply the exact quantization to study the ro-vibrational energy states 
of some diatomic molecules. To begin, we write Schrodinger equation with The Deng-Fan 
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diatomic molecular potential as 


d^R 


neir 


dr^ 


2fi 

¥ 




D 1 


2/ir2 


Rni{r) = 0, 


(5) 


where n and I denote the radial and orbital angular momentum quantum numbers, r is 
the internuclear separation of the diatomic molecules, and E^i is the bound-state energy 
eigenvalues. The /i and V (r) represent the reduced mass and interaction potential respec¬ 
tively. It is well-known that the equation of form ([5]) is an exactly solvable problem for 
s-wave. But for wave states, the problem is not analytically solvable. To obtain the 
bound state solutions, we must therefore resort to using an approximation similar to other 
works 48l-l5l| to deal with the centrifugal term or alternatively, to solve numerically. It is 


noted that, for short potential range, the following formula is a good approximation to the 
centrifugal term: 


a 


do + 


(gar _ X)2J ’ 


( 6 ) 


where dg^^. It should be noted that this approximation reduce to the one used by Dong 


and Gu [29|, Dong et al |23[, and Qiang and Dong [24 when do = 0. Now if we consider 
the approximation ([6]) and equation ([5]), we can hnd the effective potential as 


Veffig) = D 


' - 1 - l)a^doh’^ 

2fi 




2/i 


2/i 


(7) 


where we have introduced a new transformation of the form g = ■ Now, let us begin 

the application of the quantization rule to study the potential 14//(p). To perform this task, 
we have to hrst calculate the turning points Qa and Qb determined by solving V{g) = Eni- 
Thus we have 


Qa — 


R. 

2R 


- 4i?(P - Eni) 


2R 


Qb — 


Q ^jQ^-m{P- E, 
'2R^ 2R 


ni 


with k{g) = 


( 8 ) 


where k{g) is the momentum between the two turning points ga and gb- We have also 
introduced three parameters P, Q and R for mathematical simplicity. These parameters 
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are as follows: 


P = D + 


+ l)a‘^doh^ ^ £(£ + l)aV ^ ^ t{i + l)a^don^ ^^2 


2yU 


-, Q = 


2yU 


2Dh and R = 


2yU 


+ Db\ (9) 


The non linear Riccati eqnation for the gronnd state is written in terms of the new variable 
g as 


—a^(l + g)(pQ{g) + <Po{g) + [-E'o£ — Veffi^)] 0o(^) — 0. (10) 

Thus, since the logarithmic derivative 0o(^) for the ground state has one zero and no pole, 
therefore we assume the following solution for the ground states 


00 (p) = A + Bg. 


( 11 ) 


On substituting equation flTTl) into ffTOj) and then solve the non-linear Riccati equation, we 
obtain the ground state energy as 

= ( 12 ) 

Also, A and B are as follows: 




jjL Q — R B 
¥ B 


and 


B 



(13) 


It is worth to be noted that we only choose the positive sign in front of the square root 
for B. This is because the logarithmic derivative 0o(p) will decrease exponentially, which 
is required physically. Let us now calculate the quantum correction. For this purpose, we 
utilize the integrals given by appendix A, and we obtain 



K{r)^^dr 


"Qb 


fco(g) 0o(g) ^^_l 

Iga «P(1 + Q) 0o(p) « 


2pi? 




[A + B]g[g-{^)]dg 
Bg{l + g)\/{ga - g){g- gt) 


2fiR 


"Qb 


dg 




' Qa 


2/iR 




1 + 


\/{Qa - g){g- gb) 

Bh 


(f-i)(^ + i) r(^) 

1 + p g 


+ 1 


(14) 
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Furthermore, the integral of the momentum k{r) can be found as follows: 



r*’ k{g) , _ 1 /2/ui? x/(pa - p)(p- Qb) 

4 ag{l + g) ^ g{l + g) 


TT 

a 

71 

a 


2fiR 


r 


— 1 — gagb 


R — Q -\- P — Eni 
R 


P-E. 


ni 


R 


(15) 


where we have used an appropriate standard integral in the appendix A. Now by combining 
the results obtained by equations flT^ and flTSl) with equation equation ([3]), i.e. 


TT 

a 


R-Q + P-E, 


ni 


R 


- 1 - 


P-E. 




R 


Ntt H— 
a 



\ . VWR 

Eh 


(16) 


the energy eigenvalues spectrum can then be found as 


Erie — D{h + 1)^ + 


E{E+l)a^n^do 

2fi 


a? {B + an) 2nDb{b + 2) 
2p _ 2a 2h?a {B + an) 


(17) 


Let us now obtain the corresponding wave function for this system. For this purpose, we 
introduced a new transformation of the form z = G (e",0) in equation ([5]) which 
maintained the hniteness of the transformed wave functions on the boundary conditions 
we have 


_^ 2 (RRni{z) , JRneiz) , fU + Vz + Wz^ 
z ; A z ' ~t” 


dz^ 


dz 


{i-zY^ 


Rni{z) = 0 


2/i 


with W = -^ {Eni -D)- i{i + l)do, 


1 / = 


a'^h 
A^bD Afi 


a 


2^2 am^ 


{Eni-D)-e{e+l){l + 2do) 


and U — ^ {Eni ~ D) -^ 2 + 2) + l{i. + l)(io- 

h a^n 


(18) 


Furthermore, equation ffTSj) is transformed into a more convenient second-order homoge¬ 
neous linear differential equation via the following transformation 


Rni{z) = ZP{1 


-‘G. 


nRZ) 


with 


P=2+'^ 


■w. 


(19) 
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Substituting of equation (fT^ into equation ffTSD . we can find 


G"ni{,z) + 


(2p+i)-q2p+| + i) 

41 - z) 


+ Gne{z) 


41 - 4 


= 0. 


( 20 ) 


The solution to the above second order differential equation can be expressed in terms of 
the hypergeometric function as: 


G„Az)= oFGC. 




The following notations: 

C = P + — - V-U, p = p + — + y/-U 7 = 2p + 1 (22) 

a a 

has been introduced in equation (1^ so as to avoid mathematical complexity. By con¬ 
sidering the finiteness of the solutions, Gni{z) approaches infinity unless C is a negative 
integer. This is an indication that Gni{z) will not be hnite everywhere unless we take 
P + ^-V^ = —n. Thus, p given by equation fl22|) can be re-written in terms of this 
condition, and hnally we can write the wave function as 

H ^ ;2p + 1; z). (23) 

a J 

Nni is the normalization constant. 


Rni{z) = Nniz'P{l - z)o‘Gni{z) = z^{1 - z) 2Fi{-n,2 


IV. RESULTS AND CONCLUSION 


By using the known spectroscopic values in Table |T] we obtain the energy states of few 
selected diatomic molecules for various vibrational n and rotational ^ angular momentum 
as shown in Tables Hand [ml These diatomic molecules are H 2 {X^T.+), C'0(Xm+), 
i7F(X^E+), 02(X^E+) and Oj (X^IIg). The spectroscopic parameters are taken from 


the work of Kune-Gordillo-Vazquez ^ and Oyewumi et al. We also applied following 
conversion = /i x 6.0221415 x 931.494028e6/c^ with he = 1973.29eUA throughout 

our numerical computation. 


In order to test the accuracy of our results,we give a numerical comparison of our ob- 
















TABLE I: Model parameters of the diatomic molecules studied in the present work. 


Molecules(states) /^/lO ^^{g) ’’e(A) D{cm a(A) 


H 2 (xis+) 

0.084 

0.741 

38318 

1.9506 

CO (Vis+) 

1.146 

1.128 

90531 

2.2994 

HF (XiE+) 

0.160 

0.917 

49382 

2.2266 

0 

CO 

1.337 

1.207 

42041 

2.6636 

0 + (xm,) 

1.337 

1.116 

54688 

2.8151 


tained energy spectrum with the ones H 2 and CO obtained previously in 

the literature, for the shifted Deng-Fan potential (i.e. Deng-Fan potential shifted by the 
dissociation energy D). As it can been seen from Table [ITl our approximate results are 
in excellent agreement with the ones obtained previously by other authors and via other 
approach. We also proceed furthermore to obtain the ro-vibrational energy spectrum (in 
eV units) for the 5 selected diatomic molecules. 

It is worth to be noted that the advantage of the method presented in this study is that, 
for exactly solvable quantum system, it enable one to find the energy spectrum directly in a 
simple way. Finally, we recommend that the extension of this method to obtain the bound 
state solution some other potentials such as: Wood-Saxon potential, Yukawa potential, 
Hellmann potential, etc. 


Appendix A: Some useful standard integrals 


rrB I 

Jr A y/(r - rA){rB - r) 


(Al) 


j-VB I 

Jr A {a + br) 


71 

{a + brB){a + brA) 


(A2) 



-\/{xJ^VA){rB^^dr 


TT 

2 


{rA + Tb) 


TTv'rlri'. 


(A3) 
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TABLE II: Comparison of the bound-state energy eigenvalues —{Eni — D)(eV) of H 2 and CO 
molecules for various n and rotational I quantum numbers in Deng-Fan diatomic molecular po¬ 
tential. 


n I Present {Hi) AIM N-U [12] (Hi) Present (CO) AIM (CO) N-U[12j (CO) 


0 

0 

4.400978574 

4.394619779 

4.39444 

11.08220218 

11.08075178 

11.08068 


5 

4.183429818 

4.176618048 

4.17644 

11.07491243 

11.07253985 

11.07247 


10 

3.629702678 

3.621838424 

3.62165 

11.05547579 

11.05064581 

11.05057 

5 

0 

1.764237205 

1.758451567 

1.75835 

9.712515712 

9.688146187 

9.68809 


5 

1.623366051 

1.617410615 

1.61731 

9.705448303 

9.680226284 

9.68017 


10 

1.266759570 

1.260451640 

1.26034 

9.686604517 

9.659110919 

9.65905 

7 

0 

1.082609838 

1.077636993 

1.07756 

9.191396379 

9.159164003 

9.15911 


5 

0.966873845 

0.961814782 

0.96174 

9.184417149 

9.151359661 

9.15131 


10 

0.675048455 

0.669844065 

0.66976 

9.165808494 

9.130552425 

9.13050 
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TABLE III: Ebergy spectra of Deng-Fan diatomic molecular potential for 
CO HF O 2 and O^ (JA^IIg) molecules for various n and rota¬ 

tional i quantum numbers. 


n i 

H 2 (xis+) 

CO (XiE+) 

HF (XiE+) 

O 2 (x3e+) 

0 + {X^Tlg) 

0 0 

0.365141571630 

0.14496907267 

0.29664475468 

0.10204385069 

0.12341159004 

10 

0.998213655071 

0.43017899550 

0.84879137870 

0.30157505570 

0.36506351446 

1 

1.011736491223 

0.43071963141 

0.85396890484 

0.30203765756 

0.36559505801 

2 0 

1.582601950199 

0.71191080079 

1.37474141491 

0.49742226304 

0.60257089273 

1 

1.595035333506 

0.71244756155 

1.37972360389 

0.49788160541 

0.60309875940 

2 

1.619836697120 

0.71352107116 

1.38968398285 

0.49880028219 

0.60415448370 

3 0 

2.107302268265 

0.98963235238 

1.86975698501 

0.68912972889 

0.83540994827 

1 

2.118708590334 

0.99016524991 

1.87454776989 

0.68958581969 

0.83593414705 

2 

2.141459973433 

0.99123103314 

1.88412543561 

0.69049799334 

0.83698253554 

3 

2.175434452445 

0.99282967815 

1.89848217862 

0.69186623393 

0.83855509594 

4 0 

2.575691305740 

1.26335597760 

2.33445776235 

0.87670753899 

1.06359208920 

1 

2.586127919535 

1.26388502379 

2.33906095729 

0.87716038610 

1.06411262902 

2 

2.606943719809 

1.26494310432 

2.34826353550 

0.87806607249 

1.06515369967 

3 

2.638024369517 

1.26653019550 

2.36205787817 

0.87942458223 

1.06671528329 

4 

2.679199648264 

1.26864626172 

2.38043256798 

0.88123589159 

1.06879735299 

5 0 

2.990875171736 

1.53309394464 

2.76944475234 

1.06016573964 

1.28712867897 

1 

3.000394872082 

1.53361915128 

2.77386405693 

1.06061535100 

1.28764556873 

2 

3.019380399287 

1.53466955280 

2.78269894439 

1.06151456584 

1.28867933936 

3 

3.047724487910 

1.53624512555 

2.79594197541 

1.06286336840 

1.29022997301 

4 

3.085267435478 

1.53834583407 

2.81358200193 

1.06466173494 

1.29229744290 

5 

3.131798515096 

1.54097163107 

2.83560417964 

1.06690963399 

1.29488171325 

6 0 

3.355714996849 

1.79885846293 

3.17530096333 

1.23951433788 

1.50603103656 

1 

3.364366448329 

1.79937984182 

3.17953996677 

1.23996072135 

1.50654428521 

2 

3.381618777060 

1.80042258783 

3.18801433936 

1.24085348039 

1.50757077359 

3 

3.407371283026 

1.80198667737 

3.20071681651 

1.24219259927 

1.50911048386 

4 

3.441474029835 

1.80407207520 

3.21763651189 

1.24397805443 

1.51116338944 

5 

3.483729160505 

1.80667873414 

3.23875892977 

1.24620981439 

1.51372945463 

6 

3.533892617999 

1.80980659543 

3.26406598158 

1.24888783984 

1.51680863506 

7 0 

3.672849758377 

2.06066168353 

3.55259204944 

1.41476330171 

1.72031043670 

1 

3.680677877082 

2.06117924630 

3.55665423437 

1.41520646508 

1.72082005302 

2 

3.696286607610 

2.06221436015 

3.56477505504 

1.41609278397 

1.72183927685 

3 

3.719581354261 

2.06376700163 

3.57694741679 

1.41742224263 

1.72336809040 

4 

3.750421261686 

2.06583713552 

3.59316068801 

1.41919481768 

1.72540646704 

5 

3.788620441740 

2.06842471492 

3.61340071207 

1.42141047769 

1.72795437139 

6 

3.833949577906 

2.07152968119 

3.63764982344 

1.42406918344 

1.73101175897 

7 

3.886137881655 

2.07515196399 

3.66588686789 

1.42717088787 

1.73457857662 
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